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ABSTRACT 


As part of an on-line simulation analysis system being developed 
at the Naval Postgraduate School there was a need for a GPSS-like 
simulator which could be run under control of a FORTRAN program. 
The objective of the work done for this thesis was to design an informa- 
Hen structure appropriate for such . simulator and to write the sim- 
ulation routine. This thesis gives a detailed description of the 
information structure, telling what it must look like when the routine 
is called and what it contains when the routine is through. Another 
section dis cusses the procedure by which the routine performs the 


simulation. Program listings anda complete sample problem are 


included, also. 
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I. INTRODUCTION 


A research project at the Naval Postgraduate School has as its 
goal the development of a system with which simulation analyses can 
be performed through natural language dialogue with a computer [ 1}. 
Currently, the capability exists for translating an English description 
of a queueing problem into an Internal Problem Description (IPD) and 
for translating an IPD into an equivalent English description and into 
a GPSS simulation program [1, 2, 3]. ‘The simulation can then be per- 
formed by running the GPSS program [4,5]. 

It would be desirable a the system er developed to provide 
for analyst-computer interaction during the simulation, but the current 
mode of operation does not allow this, because of the GPSS program 
having te be run separately. In order to make this interactive 
capability possible, it was decided to ieaplerneet a GPSS-like simulator 


that could be run under control of the current FORTRAN program. 


A. OBJECTIVE 

The objective of the work done for this thesis, then, was to 
design an information structure appropriate for the simulator and 
to write a simulation éedtine. A design goal was to make the informa- 
tion structure as compact and self-contained as possible. This was 


accomplished by using one singly-subscripted variable common to 








the simulation routine and the calling program to hold all of the 
information about the model. 

The information in this array is structured in basically the same 
manner as the "internal tables'' of the GPSS simulator, and the 
algorithm of the simulation routine is basically the same as that of 
GPSS [5]. Although this routine is intended for use with the system 


described earlier, it could be called from any FORTRAN program. 


B. ORDER OF REPORTING 

Section II of this thesis describes the information structure in 
detail. It tells what must be there when the routine is called and what 
can be found there when it is through. Section III discusses the 
procedure by which the routine performs the simulation. Section IV 
gives conclusions and recommendations. A complete sample problem 


is contained in Appendix A. 





Il. THE INFORMATION STRUCTURE 


In the Simulation Routine all information about a model is held in 
a one-dimensional array called the X-Vector, each element of which 
is a four-byte word in the IBM 360. The allocation of elements in 
this vector is completely flexible, and, except for the first 32 elements, 
depends entirely on the particular model being specified. For exam- 
ple, if only one Storage is.used in the model, then only enough 
elements for one Storage entity need be allocated, and these elements 
can be located almost anywhere inthe array. The first 11 elements 
of the X- Vector are used for certain attributes that are needed for all 
models, and the next 21 elements furnish information about what is in 
the remainder of the array. Currently, this array is DIMENSIONed 
to have 500 elements, but this can be altered by changing 2 numbers 
in the FORTRAN program. | 

With the exception of Transactions, the X-Vector elements for 
all entities must be assigned and initialized prior to the eating of 
the Simulation Routine. Each type of entity (except Transaction and 
Savevalue) has 2 biestices in the first 32 elements of the X-Vector to 
specify the quantity of that type of entity in the model and to furnish 


a@ pointer to the directory for that type of entity. The directory, then, 


contains a pointer to each entity of that type. Each of these pointers 








actually points to the element immediately preceding that entity. The 
locations of the directories and the entities are arbitrary, except that 
they must be after element 32 and before the Transactions. 
Savevalues do not require a directory, since each one uses only 
one X-Vector element, and these elements must be contiguous. For 
transactions all that need be specified is a pointer (X25) to the begin- 
ning of the area in the X-Vector where the Transactions are to be. 
The Simulation Routine initially allocates elements for each Trans- 
action and chains these saianied to form a list of unused Transactions. 
In this section of the thesis the X-Vector layout for each type of 
entity is described in detail. Also, the information required in the 
first 32 elements is specified. It is stated what must be in the X- 
Vector prior to calling the Simulation Routine and what will be found 


there after the routine is executed. The sample problem in Appendix 


A furnishes examples of most of the points discussed in this section. 


A. FIRST 32 ELEMENTS OF THE X-VECTOR 

These elements contain attributes and pointers which ae major 
role in controlling the simulation. The following subsections describe 
the manner in which these elements must be initialized by the calling 
program. Where a value is changed during the course of the simula- 
tion, it is so indicated. If an entity is not defined in a model, then 


that entity's associated element need not be initialized to any value, 


since the Simulation Routine will never reference it. 








1. Mode (X1) 
This element is initialized to the values 1,2,3, or 4, 


according to the following table. 








MODE GPSS 
i Simulate /Start 
2 Reset/Start 
3 Clear/Start 
4 Start 


« 


This element performs the same function as the Start, Reset/Start 
and Clear/Start cards in GPSS. If the Mode is initialized to 1, the 
Simulation Routine computes slopes, sets the clocks to zero, creates 
a chain of unused Transactions, and places one Transaction on the 
FEC for each GENERATE block a the model. If the Mode is 
initialized to 2 or 3, then the model entities are reset or cleared the 
same way as by the Reset card or Clear card in GPSS. The simula- 
tion then continues from where it left off when it returned control to 
the calling program. If the Mode is initialized to 4, then the calling 
program must set the Termination Count to some new value, and 
return control to the Simulation Routine where the simulation will 
continue from where it left off. 
2. Termination Count (x2) 
The Termination Count serves the same purpose as the first 


argument of a GPSS Start card. It must be initialized to some integer 


value each time the Simulation Routine is called. Then during the 








ae 


simulation, whenever, a Transaction enters a TERMINATE block, this 
value is decremented by the value of argument A of the TERMINATE 
block. When the Termination Count reaches zero, the simulation 
stops and control is returned to the calling program. 
3. . Random Number Seeds (X3-X10) 
These elements are initialized to the same values used ona 
GPSS RMULT card. The values will change during the course of the 
simulation, since the random number generator stores the current 
random number for each pemnegen in these elements. 
4. Clock Time (X11) 
This element corresponds to the GPSS relative clock. 
5. Number of Storages (X12) 
This element is initialized to the number of Storage entities 
sieli in the model. 
6. Storage Directory Poses (X13) 
This element is initialized to the subscript of the element 
immediately preceding the Storage directory. 
7. Number of Queues (X14) 
This element is initialized to the number of Queue entities 
used in the model. 
8. Queue Directory Pointer (X15) 


This element is initialized to the subscript of the element 


immediately preceding the Queue directory. 








9. Number of Functions (X16) 


This element is initialized to the number of Function entities 


defined in the model. 
10. Function Directory Pointer (X17) 
This element is initialized to the subscript of the eect 
immediately preceding the Function directory. 
11. Number of Blocks (X18) 
This element is initialized to the number of Block entities 
defined in the model. ; 
12. Block Directory Pointer (X19) 
This element is initialized to the subscript of the element 
immediately preceding the Block directory. 
13. Number of Tables (X20) 
This element is initialized to the number of Table entities 
defined in the model. 
14. Table Directory Pointer (X21) 
This element is initialized to the subscript of the element 


immediately preceding the Table directory. 


15. Number of Savevalues (X22) 


This element is initialized to the number of Savevalue entities 


used by the model. 


16. Savevalue Pointer (X23) 


This element is initialized to the subscript of the element 


immediately preceding those to be used for Savevalues. ~ 
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17. Rig wbee ae —— (X24) 

This value is initialized to the number of parameters each 

Transaction in the model must have. 
18. Transaction Pointer (X25) 

This value is initialized to the sdiisixipt of the Sosa 
immediately preceding the Transaction entities. After all other 
entities have been assigned X-Vector locations, this element would 
normally be given the highest subscript assigned to those entities. 
The program will then sania all the remaining X-Vector elements to 
the Transaction entities. 

19. Pointer to First Transaction on Current Events Chain (X26) 

This element is initially set equal to zero. Thereafter, it 
will contain the pointer to the Gest Transaction on the Current Rvanes 
Chain (CEC), or zero is there is no Thensacean on the CEC. _ 

20. Pointer to Last Transaction on Current Events Chain (X27) 

This value is initially set equal to zero. Thereafter, it 
contains a pointer to the last Transaction on the CEC, or serast there 
is none. 

21. Pointer to First Transaction on Future Events Chain (X28) 


This value is initially set equal to zero. Thereafter, it 


contains a pointer to the first Transaction on the Future Events 


Chain (FEC), or zero if there is none.: 





22. Pointer to Last Transaction on Future Events Chain (X29) 


This value is initially set equalto zero. Thereafter, it will 





contain a pointer to the last Transaction on the FEC, or zero if there 
is none. : 
23. Error Code (X30) 
This element is initially set equal to zero. If during the 
7 simulation, an error is detected by the routine, an error code will 
be placed in this element and control will be returned to the calling 
program. “These codes chansapend to the error numbers in GPSS, 
and are listed in Figure 2. 
24. Number of Variables (X31) 
This element is initialized to the number of Variable entities 
defined in the model. | 
25. Variable Directory Pointer (X32) 


This element is initialized to the subscript of the element 


immediately preceding the Variable directory. 


B. STORAGES 

Each Storage entity defined in the model must be assigned 10 
contiguous 6 Wewter elements as illustrated in Figure 3. These values 
have the same meaning as in.GPSS., The calling program must 
initialize the second element for each Storage entity to the Storage's 
capacity and all other elements, meee the IPD pointer! to zero. 


~ 


1 currently, IPD pointers are completely ignored by the Simulation 
Routine and can simply be initialized to zero. 
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The Simulation Routine will make entries in the elements of a Storage 
each time a Transaction passes through an ENTER or LEAVE block 
that refers to that Storage. If the amount of Storage units required 
by a Transaction is less than or equal to the amount available, then 
the current contents is incremented by that amount and the nuntlge 
of units available is decremented by that amount. 

The cumulative time integral is an 8-byte floating point sum of 
the number of clock units that the Storage was occupied, weighted by 
the number of units in the Seowape during each interval that occupancy 
was constant. The integral is updated by the Simulation Routine each 
time the Storage contents changes. The.cumulative time integral can 
be used to compute Storage utilization, which is the cumulative time 
integral divided by the product of she clock time and the Storage 
capacity. | 

The last clock time the Storage changed status contains the clock 
time at which the time integral was last updated. The clock time is 
placed in this element each time the Storage changes status, and 
when the Simulation Routine is called using mode 2 (Reset) or mode 3 
(Clear). | 

Entry count is incremented every time a Transaction enters an 
ENTER block. The fon poitude of the increment is specified by 
argument B of the ENTER block. 


Maximum contents is set equal to current contents each time a 


Transaction enters an ENTER block and the new current contents is 


greater than the previous maximum contents. 
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The remaining 3 elements contain 5 pushdown delay chain 
pointers and the IPD pointer. These delay chain pointers are 
associated with the following Storage conditions, (1) full, (2) not 
empty, (3) empty, (4) not full and (5) reduction in contents. Each 
time a Transaction enters an ENTER block or a GATE block and 
fails to advance because the Storage condition associated with that 
block is not met, the Transaction is placed in the appropriate delay 
chain and the delay chain pointer associated with that condition is 


set to point to the Transaction. 


C. QUEUES 

Each Queue entity defined in the model must be assigned 8 
contiguous X.Vector elements as illustrated in Figure 4. These 
values have the same meaning as in GPSS. The calling program 
must initialize all Queue elements, except the IPD pointer, to zero. 
All future entries in these elements will be made by the Simulation 
Routine. 

The last clock time the Queue changed status contains he tint 
clock time at which the cumulative time integral was computed. The 
clock time is placed in this element each time the Queue changes 
status and when the Simulation Routine is called in mode 2 (Reset) or 
mode 3 (Clear). 

Total entries is incremented each time a Transaction enters a_ 
QUEUE block which refers to the Queue. The magnitude of the 


increment is specified by argument B of the QUEUE block. 


a7 
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Number of zero-delay entries is incremented whenever a trans- 
action enters and leaves the Queve with no time delay. The magnitude 
of the increment is given by the value of argument B of the QUEUE 
block. 


The other elements are similar to those of the Storage entities. 


'D. BLOCKS 


Block entities defined in a model describe the various actions 
which are to take place during the simulation. Each Block may have 
associated with it from 0 to 7 arguments. In general, these Blocks 
have the same meaning as Blocks in GPSS, and the arguments perform 
the same function also. 

Figure 5 is a general layout of the elements of a Block entity. 
From Z to 9 contiguous X-Vector elements are required for each block, 
depending on the number of arguments specified. The calling program 
must initialize these elements prior to calling the Simulation Routine. 
Throughout the simulation, these values will remain unchanged. 
Figure 6 is a list of Block type numbers, Block names and sheaciiend 
arguments for each Block. Each argument may be a positive integer 
constant or any other Standard Numerical Attribute (SNA). A list of 
SNA's and their associated numerical codes is provided in Figure 7. 
SNA's have the same meaning here as they do in GPSS. 

Each X-Vector element assigned to a Block entity contains two 


pieces of information, with the exception of an argument which is an 


ail | 








integer constant. Each piece of information is represented as a 
2-byte integer value. The lst element of a Block definition contains 
the Block type number in the first half of the element and the IPD 
pointer in the second half. The first half of the element ofa 
GENERATE block would have a 1 init, an ADVANCE iackeak Zs 

an ASSIGN block a 4and soon. This value is used as the index 
number in a "Computed GOTO" statement to execute the routine 
associated with the block type and corresponds to the ''Execution 
Routine Mask" in GPSS. : 

If there are 1 or more arguments associated with the Block 
being defined, then the next element will contain argument A, the 
following element will contain argument B, and so forth. If the 
argument is a positive integer ——" then the whole element will 
sorttnde the constant. All other SNA's are represented in two parts. 
The first half of the element contains the negative of the SNA name 
code listed in Figure 7. The negative value is used to distinguish 
between constants and other SNA's. The SNA name ee used as 
an index ina ''Computed GOTO" statement to execute the evaluation 
routine associated with that SNA name. The second half of the element — 
is the entity number, if it is positive. If the second half is negative, 
indirect addressing is indincaea i.e., the absolute value of the second 


half of the element is the parameter number of the current Transaction 


which contains the entity number of the SNA. The following examples 











show GPSS Blocks and arguments and the equivalent coding required 


for the Simulation Routine. 


GPSS X.Vector Block Representation 





GENERATE V1,2 


ADVANCE 6, V¥4 





The first part of the last element of a Block definition must be a 
-33 to signify to the Simulation Routine that the remaining arguments 
for this Block are zero. The second half of the last element is either 
zero or contains a numerical mode code which is used to define the 
er operation required for a TEST,GATE, OR SELECT pbk. In 
GPSS this information is givenas a modifier on the Block name. A 
table of mode codes appears in Figure 8. The Block routine uses 
this code as an index to a ''Computer GOTO" statement within the 
routine associated with the Block type. The following examples 


illustrate the use of the mode code. 


GPSS X-Vector Block Representation 





TEST®5,Q1,6 
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X-~Vector Block Representation 








caTE Gn) 2, P*3 
SELECT (an) A, 2, Sy, 95940 





E. TABLES 

Each Table entity requires a basic X-Vector element assignment 
of 9 contiguous elements. In addition, one element is required for 
each frequency class. All elements correspond to the Table allocation 
in.GPSS. Figure 9 is a general layout of a Table definition in the 
X-Vector. One major difference between the Simulation Routine and 
GPSS is that in the X-Vector the SNA to be tabulated appears as 
argument A of a TABULATE block instead of being included in the 
Table definition. The calling program must initialize element 6 to 
the width of the frequency classes, element 8 to the number of 
frequency classes, element 9 to the upper limit of the lowest interval 
and all the other elements to zero. 

The Simulation Routine will modify the elements each time a 
Transaction enters a TABULATE block which references the TABLE 


in Argument B. Sum of Arguments entered in Table is a signed 
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floating point value that is incremented by each argument value that 
is entered inthe Table. The sum is used in the computation of the 
mean argument. Sum of squared values of the arguments entered in 
the Table is a floating point value that is incremented by the square 
of each ee value that is entered inthe Table. The ites is 

used in the computation of the standard deviation of the Table argument. 
Sum of weighted value of arguments is a floating point value that is 
incremented by the product of the argument and the number of units 
specified by the Cegtmaeead weighting factor of the TABULATE block. 
The sum is used in the computation of a weighted mean. Sum of 
weighted squares is a floating point value incremented by the product 
of the square of the argument and the number of units specified by the 
argument C weighting factor of the TABULATE block. The sum is 
used in the computation of a standard deviation of the weighted 
argument. 

Number of entries in Table is incremented each time a Trans- 
action enters a TABULATE block referencing the Table. ‘The 
magnitude of the increment is specified by argument C of the 
TABULATE block. Sum of overflow arguments is a Fisting point 
value that is incremented by each argument value which lies in the 
last frequency interval of tte Table. This sum can be divided by 
the number of entities in the last interval to obtain "Average Value of 


Overflow". A Frequency Interval corresponding to the argument is 


incremented by the value of argument C of the TABULATE block. 








F, FUNCTIONS 

Each Function entity defined in the model must be assigned 6 
contiguous elements for basic information, plus additional contiguous 
elements for independent variable (xy values ifany, Function (Y) 
values, and slope values, ifany. Figure 10 is a general layout ofa 
Function as it would appear in the X-Vector. 

The Argument A SNA must be initialized by the calling program 
to the SNA that will be used as the search value. It can be any SNA 
allowed by the Simulation Routine, as listed in Figure 7. The next 
element must be initialized to a -23 in the first half and the IPD 
pointer in the second half. The third Kuaaider should be initialized to 
zero. The first half of element 4 must be initialized toa -25. These 
last two elements are used by the simulation routine to evaluate the 
function values when they are not simply constants (i.e., when the 
mode is E or M). : 

Argument B, Part 1 (Mode) and Part 2 (Number of Points) must 
be initialized by the calling program in the fifth element. These 
values have the same meaning as in GPSS. The mode is a numeric 


code which is related to the GPSS mode as indicated in the following: 


GPSS Argument B-Part 1 (Mode) 








tg 


Part 2 of this element is the number of points which define the 
Function. 

The next element (the sixth) is initialized by the calling program. 
Part 1 (X value Integer/Decimal Indicator) must be set equal to 0 if 
the X values are integer and 1 if the X values are floating — 

Part 2 (Y value Integer/Decimal Indicator) must be set equal to 0 if 
the Y values are integer and 1 if the Y values are floating point. 

The X values and Y values must be initialized by the calling 
program to the values abit define the Function. If the Function 
mode is continuous (Mode 1) then the calling program must reserve 
elements immediately following the Y values for slopes. The number 
of elements required is equal to the number of Function points. If 
the Function mode is discrete numerical or discrete attribute (Mode 2 
we 3), then no elements are required for the slopes. If the Funetten 
is a list numerical or list attribute (Mode 4 or 5), then no elements are 
required for slopes or X values, since no:s aairch is performed. For 
these list mode Functions, the Y values must begin immediately in 
the seventh element. In those cases where slopes are needed they 


are computed in floating point by the Simulation Routine. The following 


examples illustrate the X-Vector representation of various GPSS 


Function definitions. 











GPSS 





FUNCTION RNI, C2 
, 0.3,4.0/1.0,11.0 


FUNCTION P3, D2 
6, 4/8, 10 


FUNCTION V*4, M2 
1, Q#1/2, Q¥2 


X-.Vector Function Representation 
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GPSS _ X-Vector Function Representation 





FUNCTION RN3, E3 
0.0, V*1/. 542, 15739688/ 
1.0, V10 





G. VARIABLES 

Each Variable defined in the model must be completely specified 
by the calling program. A Variable definition requires contiguous 
X-Vector elements, with the number of elements depending on the 
complexity of the Variable. Each Variable is defined in ''Polish 
spent wine notation. The operand elements may be any SNA, including 
an integer constant, in the format discussed earlier. The operator 
elements must have a -30 in the first half and an operation code in 
the second half. This code is used to indicate which arithmetic 
operation is to be performed on the previous two evaluated arguments. 


The arithmetic codes are listed below. 


Code Arithmetic Operation 


Add 

Subtract 
Multiply 
Divide 

Divide Modulo 


Ob WN Re 








The last element of a Variable definition must be initialized to a 
-25 in the first half. The following example shows two different ways 
in which the same Variable could be defined. 
GPSS X-Vector Variable Definitions 


VARIABLE 16+4* FN2 





H. TRANSACTIONS 

Each Transaction entity requires a number of contiguous 
X.Vector elements which is equal to 8 plus the number of parameters 
assigned each Transaction (according to the value of X24). Figure 11 
shows the layout of a Transaction entity. The calling program need 
aay initialize Transactions in any way; the Simulation Routine will 
make all entries in these elements during the simulation. 

When a Transaction is inthe FEC, its first element contains the 
location of the previous Transaction in the FEC in the first half and 
the location of the next Transaction in the FEC in the second half. 
When a Transaction is in the CEC, the first half of this element con- 
tains the delay chain address of the next Transaction in the Chain, 
and the second half contains the current Block number. The second 
element contains similar values when the Transaction is in the FEC, 


except that there is no Transaction which is in the FEC and in a delay 


chain at the same time. The third element contains the clock time that 





the Sacrweetiion is scheduled to leave the FEC. Each time the Trans- 
action enters an ADVANCE block, this element will be set to some 
future time if the time delay is positive. The fourth element is the 
Absolute Clock time which is assigned each time the Transaction 
passes through a MARK block. | 

The first half of the fifth element contains the number of the 
next Block that the Transaction is ed pcicien to enter and the second 
half contains the priority. The sixth element has two logical indicators 
in the first half anda — number, when the Transaction is ina 
Queue, inthe second half. The first indicator is not currently used 
and the second,is in the scan status oe which is used by the 
Scan routine to determine if the Transaction is eligible for attempting 
to move it into its next Block. The next element is the time the 
tid Seae entered the Queue whose number is in element 6: The 
eighth element contains the IPD pointer in the first half and the 
Transaction number in the second half. Transactions are numbered 
consecutively as they are created in GENERATE blocks. The- 


remaining elements are the Transaction parameters, which are used 


in the same manner as in GPSS. 


¢ 


I. ASAMPLE PROBLEM 
A sample problem is presented in Appendix A to show a complete 


X-Vector, both before and after execution of the routine. In the 


Appendix an equivalent GPSS program is given first. Then there is a 





listing of the cards used to initialize the X-Vector, followed by a print- 
out of the X-Vector, produced by the Simulation Routine after papeied 
Finally there is a listing of the short FORTRAN program used to 
initialize the X-Vector and call the Simulation Routine. 

The input listing can be seen to consist of four columns of numbers. 
Column 1 is the X-Vector subscript of the element to be loaded. 
Column 2 is 0 if the information is 4-byte integer, 1 if the information 
is two 2-byte — and 2 if the inforrhation is a floating-point 
value. When column 2 ccieSe 1, column 3 contains the first 2-byte 
integer. When column 2 equals 0 or 2, column 4 contains the 4-byte 
integer or floating-point value. 

The X-Vector print-out can also be seen to consist of four columns 
of numbers. Column 1 is the X-Vector subscript of the element, col- 
umn 2 is its hexadecimal representation, and columns 3 and 4 ave 


the two 2-byte integer halves of the element. When the element has 


a floating-point value, it is printed in decimal form with an exponent. 
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Im. THE PROCEDURE 


The Simulation Routine is a FORTRAN subprogram which 
performs a simulation by manipulating information in the X-Vector. 
The algorithm which it follows is essentially the same as that of 
GPSS, with the X-Vector serving the same role as the "internal 
tables'' of GPSS. The Simulation Routine consists logically of 
several routines, but physically it is just one subroutine (except for 
the Error handling routine). 

In this section of the thesis the overall flow of the Simulation 
Routine will be described first, followed by a fairly detailed discussion 
of each of its parts: Scan and Move Transaction, Argument Evaluation, 
Block Routines, and Update Clock. Finally, there are discussions of 
tracing and error handling. A listing of the FORTRAN program is 


included at the end of the thesis, with a dictionary of most of the 


program variables appearing in Appendix B. 


A. OVERALL FLOW OF ROUTINE 

When the Simulation Routine is given control, it will perform the 
appropriate initialization of model entities according to the mode 
specified in X(1). Control ie tities given to the Scan routine which will 
attempt to move each Transaction on the CEC through as many Blocks 


as possible. Each time a Transaction attempts to enter a Block, the 


Block arguments will be evaluated and the Block routine performed. 





When no Transaction can move into some next Block, control is given 
to the Update Clock routine which will increment the relative clock to 
the next event time in the FEC. All Transactions with this new clock 
time that are in the FEC will then be moved to the CEC. Control 

will again be given to the Scan routine. This cycle which is ‘Mustrated 
by Figures 12 through 14, will continue until some event within the 


model halts the simulation or an error occurs. 


B. SCAN AND MOVE TRANSACTION 

When control is passed to this routine, all Transactions with 
an event time less than or equal to the relative clock (X11) will have 
been merged into the CEC. These Transactions are chained together 
in the CEC in descending order of priority, as in GPSS. The Scan 
routine will reset the status change flag (SC FLAG) to off and examine 
the first Transaction on the CEC. If the scan status indicator (SSIND) 
is true, then an attempt will be made to move the Transaction into 
its next Block. If the SSIND is false, then some blocking condition 
within the model has caused the Transaction to become ineligible for 
movement into some next block. This condition will persist until the 
blocking condition is removed and the SSIND set to true or the 
Simulation Routine is called using mode 3 (Clear/Start) in which case 
all Transactions in the model are destroyed. 


The Scan and Move Transaction routine illustrated in Figures _ 


12 and 13 will move each Transaction on the CEC one at a time 





through as many Blocks as possible. Each time a Block routine is 
executed, the SCFLAG is checked to see if there has been a status 
change, as in GPSS. If there has been a change, the SSIND for all 
appropriate Transactions is set totrue. The Scan and Move Trans- 
action routine will then continue to move the Transaction = it 
experiences a positive time delay inan ADVANCE block, enters a 
TERMINATE block or is blocked from moving into some next Block. 
If a positive time delay is encountered in an ADVANCE block, the 
Transaction is unlinked — the CEC and merged into the FEC in 
ascending order of Block departure time (BDT). Next the SCFLAG 
is checked. If it is on, the scan starts over with the first Transaction 
in the CEC. Otherwise the next sequential Transaction in the CEC 
is processed. When no Transaction can move into a next Block, 
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control is transferred to the Update Clock routine. 


Cc. EVALUATE ARGUMENTS 

In section II it was stated that the arguments of a Block could be 
any SNA. These SNA codes may be considered to be ad eoctiies 
for a ficticious computer, indicating what routine must be performed 
to evaluate them. The following algorithm and flowchart describe in 
detail the argument evaluation procedure. 


To evaluate SNA's requires 3 stacks and three pointers as indicated 


below. 





PNTA 


INST 







current 
argument 





As illustrated in the flowchart in Figure 15 the pointers PNTA and 
PNTB are initially set to zero anna INST is set to the location (sub- 
script) of the first argument of a Block. The fullword (4 byte) value 
of the argument is then checked for sign. If it is positive the whole 
argument is considered a positive integer constant and placed on the 
OS/FOS stack as a positive integer. If itis negative, the argument 
is either an SNA listed in Figure 7 or a code (-23, -25, -30, or -33). 

All evaluated arguments will eventually appear in consecutive 
locations on the OS/FOS stack, and the cortesponding integer/decimal 
indicator (INTDEC) elements will contain a 0 or 1 indicating whether 
the value on the OS/FOS stack is an integer or floating point value, 
respectively. 

To place an evaluated argument on the OS/FOS stack, a routine 


must be performed which translates the SNA code to some value. 


The first thing which must be done is to check the second half of the 
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argument for a negative value which indicates indirect addressing. 

If the value is negative, the contents of the parameter, indicated by 
the absolute value of this number, contains the entity number portion 
of the SNA. This number is saved for later use. If the value of the 
second half is positive, it is the entity number to be used. The 
negative of the first half of the argument specifies the SNA name 
which is used as the index ina "Computed GOTO" statement. 

When the "Computed GOTO" statement is executed, the appropriate 
routine produces a value we the SNA and places it on the top of the 
OS/FOS stack. The INST pointer is then incremented to point at the 
next Block argument. This procedure continues untila -33 code is 
detected. When the first half of the argument is a -33 code, this 
signals the end of arguments for this Block. The routine associated 
swith code -33 determines whether a full seven arguments have been 
evaluated and placed on the OS/FOS stack. If not, this routine pads 
integer zeros into sequential OS/FOS and INTDEC locations after 
the last evaluated argument, up to and including the seventh location. 

1. Variable Evaluation 

The evaluation of a Variable is treated essentially as a 
recursive subroutine call upon the argument evaluation routine. The 
location of the next argument is saved on the return address stack 


(RAS), and then INST is set to the location of the first element of the 


Variable definition in the X-Vector. At this point argument evaluation 


continues with the operands of the Variable being evaluated and placed 





on the operand stack (OS/FOS). Whena Variable element is sensed 

“ with a -~30 code in the first half, the arithmetic operation indicated 
by the code in the second half is performed on the last pair of evaluated 
operands on OS/FOS. The result then replaces these top two operands 
on the operand stack. Whena Variable element is sensed a a -25 
in the first half (at the end of the Variable definition), a return is 
made by popping the top return address off RAS and putting it in INST. 
The truncated value of the Variable is left on the top of OS/FOS as 
an integer. A Variable definition may contain other Variables. The 
only restriction is that a Variable must not have itself imbedded in 
its definition at any level. 

= 2. Function Evaluation 

When a Block argument has a -22 code in its first part, the 

i irenlidtdive Routine must depart from normal argument evaluation as 
it did in the case of a Variable. First the location of the next argument 
is placed on top of the return address stack (RAS). Then INST is set 
to the location of the Ist element of the Function definition which is 
the SNA to be used as the search argument. The argument evaluation 
routine then places the value of this SNA on top of the OS/ FOS stack. 
The next argument to be evaluated will have a -23 code in the first 

. half, which will cause the mode of the function to be determined. If 


the mode indicates that the Function is not "Attribute Valued", then 


the Function value is placed directly on the OS/FOS stack and a return 


‘ 


is made to pop the location of the next argument to be evaluated off 





the RAS stack. If the Function is "Attribute Valued" then the SNA to 
be returned as the Function value must be evaluated. By placing 

the SNA into the 3rd element of the Function definition, continuation 
of the argument evaluation routine will put the value of the SNA on 
top of the OS/FOS stack. The next argument will have a 25 code 

in the first half which will cause the Simulation Routine to return and 
pop the location of the next argument off RAS. 

When locating the Function value to be returned, the same 
procedure is used for ia mode as is done in GPSS. The search 
value and the Function values can be any SNA including another 
Function or a Variable. The Variable could have as its argument 
some other Function and so on. As long as the Function does not 
have itself imbedded in the search argument or in one of the Function 


values, there is essentially no limit to the depth to which this ‘scheme 


can be carried. 


D. BLOCK ROUTINES 

After all arguments have been evaluated and their values ieee 
on the OS/ FOS stack, the Simulation Routine will execute a routine 
which corresponds to one of the Block type codes which are listed in 
Figure 6. Each routine is a section of FORTRAN code designed to 
perform the same function as the Block routines in GPSS. These 


routines have different arguments as listed in Figure 6, and each _ 


routine looks in particular elements of the OS/ FOS stack for its 





corresponding evaluated arguments. Normally, the Transaction 
being processed will proceed to the next sequential Block, except 
that certain Block routines may cause the Transaction to be delayed 


or transferred to some none sequential Block. 


E. UPDATE CLOCK 

As illustrated in Figures 12 and 14 this routine is executed when 
Transactions on the CEC can no longer move into some next Block. 
The relative clock (X11) is updated to the next most imminent event 
on the FEC, as in GPSS. All Transactions on the FEC with this 
new Block departure time (BDT) will then be moved to the CEC. 

As Transactions are moved to the CEC, they are merged into 
the chain by priority. The Transactions are chained on the CEC in 
déscending priority, i.e., highest priority number first. Each new 
Transaction becomes the last Transaction on the chain with that 
priority number. 


After all Transactions have been moved, control is returned to 


the Scan routine. 


F. TRACING VARIABLES 

The Simulation Routine has provisions for tracing the execution 
of the simulation. Through the use of the NAMELIST feature in 
FORTRAN, it is possible to set the logical trace variables and vary 


the amount of tracing output. Trace variable TR6 is common to 


both the Simulation Routine and the calling program. If TR6, which 
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has a default value of N'false", is set to true", then the amount of 
tracing output may be specified by the user. With TR6 set to ''true"’, 
the Simulation Routine will request the user to specify the logical 


value of various trace switches which control the tracing output. 


G. ERROR CODES 

Figure 2 is a list of codes for errors detectable by the Simulation 
Routine. When an error is detected, a subroutine is called. The 
arguments provided to the “subroutine are an error number from 
Figure 2, a recovery point number (presently being used to locate the 
code where the error occurred since some errors can occur in more 
than one place), and the statement number to which the error sub- 
routine is to return. The pecovery number could be used in future 


expansion of the Simulation Routine to provide a link for continuing 2 


simulation after an error has been corrected. 


i 
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Vv. CONCLUSIONS AND RECOMMENDATIONS 


A GPSS-like simulator callable froma FORTRAN program has 
been developed. This subroutine can be used by any program which 
properly initializes the X-Vector with information about the model, 
in a format similar to the "internal tables" of GPSS. Although only 
a portion of the complete GPSS capability was implemented, both the 
information structure and.the Simulation Routine are such that 
expansion is straightforward. 

It is recommended that work be done at this time to incorporate 
the simulator into the system described in Reference 1 for performing 
simulation analyses through natural language interaction with the 
computer. What is required essentially is a set of ''encoding rules" 
which will initialize the X-Vector from an Internal Problem Descrip- 
tion, call the Simulation Routine, and then report the results of the 
simulation. As the overall system is expanded to handle additional 
types of problems, the capabilities of the simulator can also be 


expanded accordingly. 


39 








TERMINATION 
COUNT 


CLOCK TIME 


NUMBER OF 
STORAGE: 


STORAGE 
DIR. PNTR. 


NUMBER OF 
QUEUES 
QUEUE 
DIR. PNTRe 


NUMBER OF 
FUNCTICNS 


FUNCTION 
Dir. ENT. 
NUMBcH OF 
BLOCKS 


BLOCK 
DIRs PNTRe 


NUMBER OF 
TABLES 


TARLS 
DIR. PNTR. 


NUMBER OF 
SAVEVALUES 


SAVEVALUE 
POINTER 


NUMBER OF 
PARAMETERS 
TRANSACTION 


POINTER 


FIRST TRANS 
CEC. PNTR. 


LAST TRANS» 
CEC. PNTR. 


FIRST TRANS. 
PEC. PNTR. 





Figure 1 X-VECTOR GENERAL LAYOUT 


o 


LAST TRANS. 
FEC. PNTR. 


ERROR CODE 
NUMBER OF 
VARIABLES 


VARIABLE 
DIR. PNTR. 


at 
ae 
Bare 
adh Sad 
ee 


STORAGE DIR. 
& ALLOCATION 


te ce 


QUEVE DIR. 
& ALLOCATION 


FUNCT. DIR. 
& ALLOCATION 


VRBL. DIR 
& ALLOCATION 


SAVEVALUE 
ALLOCATION 


es 


TRANSACTION 
ALLOCATION 


* 








BLOCK DIR. 
& ALLOCATION 


TABLE DIR. 
& ALLOCATICN 


40 








CODE MEANING 





; 2 Illegal Block type code 
45 Illegal SNA 
401 No new event in system 
425 Transaction leaving Storage by more than Storage contents 
428 Transaction leaving Queue by more than Queue contents 
433 Illegal Savevalue number 
435 Illegal Table number 
468 Number of Tenneson exceeded 
492 Illegal Transaction parameter number 
; 499 Illegal Storage number 
2 500 Illegal Queue number 
505 Negative time delay computed 
507 Illegal Function number 
509 Illegal index evaluated for list-type Function 
514 Illegal Variable number 
603 Illegal Block number 


Figure 2 ERROR CODES 2 ~ 
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Figure 3 STORAGE ENTITY LAYOUT 
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Figure 4 QUEUE ENTITY LAYOUT 
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Figure 5 BLOCK ENTITY LAYOUT 
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BLOCK TYPE CODE BLOCK NAME AND ARGUMENTS 


1 GENERATE 
A-Mean time 
B-Spread/ Function modifier 
C-Initialization interval 
D-Creation limit 
E-Priority level 





2 ADVANCE 
A-Mean time 
a B-Spread/ Function Modifier 


3 TRANSFER 
B-Next Block 


4 _* ASSIGN 
A.Parameter number 
B-SNA to be assigned 
* C.Function modifier 


: 5 MARK 
A-Parameter number 


6 ENTER ' 
A-Index number of Storage 
B-Number of units to enter Storage 


7 LEAVE 
A-Index number of Storage 
B_Number of units to leave Storage 


8 QUEUE 
A-Index number of Queue 
B-_-Number of units to enter Queue 


9 DEPART 
A-Index number of Queue 
B-Number of units to leave Queue 


: 10 TERMINATE 
A-Termination count 


: 11 * TEST 
: A-First SNA 
B-Second SNA 
C-Next Block for a false relation 


Figure 6 BLOCK TYPE CODES AND BLOCK ARGUMENTS 
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12 TABULATE 
* A-SNA to be tabulated 
* B-Index number of Table 
. * C-.Weighting factor 


13 * SAVEVALUE 
A-Savevalue number 
B-SNA to be saved 


16 * GATE 
A-Index number of Storage 
B-Next Block if false condition 


21 * SELECT 
A-Parameter number 
% B-Lower limit of entity class 
C-Upper limit of entity class 
D-Comparison value, if applicable 
E-Entity to be examined 
F.Next Block if condition not met 





*The Blocks and Arguments marked differ slightly from their GPSS 
counterparts, as discussed in Section II-D 


- Figure 6 continued 
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ANDARD NUMERICAL ATTRIBUTE (SNA) NAME CODE 


ST 


: Parameter (P) 1 
f Priority (PR) a 
Transit Time (M) . 3 
: Parameter Transit Time (MP) 4 
Storage Entry Count (SC) 5 
Storage Maximum Contents (SM) 6 
Storage Average Utilization (SR) 7 
Storage Average Contents (SA) 8 
Storage Current Contents (S) 9 
Storage Remaining Contents (R) 10 
- Storage Average Time per Transaction (ST) 11 
ES Queue Current Length (Q) 12 
; Queue Average Contents (QA) j 13 
Queue Maximum Contents (QM) : d 14 
Queue Entry Count (QC) 15 
Queue Zero Entries (QZ) - 16 
Queue Average Time per Transaction (QT) 17 
Queue Average Time per Transaction excluding Zero (QX) 18 
Table Mean (TB) 19 
* Table Entry Count (TC) 20 
Table Standard Deviation (TD) 21 
: Function (FN) | : : 22 


Figure 7 STANDARD NUMERICAL ATTRIBUTE (SNA) CODES 
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Variable (V) 
Random Number (RN) 
Clock Time (C) 


Savevalue (X) 


' Figure 7 


continued 


48 


24 


26 


27 
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CODE SELECT TEST GATE ASSIGN /SAVEVALUE 





1 E E SNE Subtraction 
2 L L SF Replacement 
3 LE LE SNF Addition 
4 GT GT SE 
5 GE GE 
6 U U 
7 MAX 
8 MIN ; 
9 SNE 
10 SF 
. 11 SNF 
12 SE 


Figure 8 MODE CODES 
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Figure 10 FUNCTION ENTITY LAYOUT 
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Figure 14 FLOWCHART OF UPDATE CLOCK ROUTINE 
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Figure 15 FLOWCHART OF ARGUMENT EVALUATION ROUTINE 
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APPENDIX A 


SAMPLE PROBLEM 


SIMULATE 
RMULT 277, 423, 715, 121, 655,531, 999, 813 
TABLE Ml,,1,2 

; TABLE Ml,,1,2 


STORAGE 1 

FUNCTION RN1, C24 

0,0/.1,.104/.2, .222/.3, _355/.4, .509/.5, -69/.6,.915/.7,1.2/ 
=75, 1. 39/. 8, 1. 6/. 84, 1. 83/. 88, 2. 12/.9,2.3/. 92,2. 52/.94,2.81/ 
95, 2.99/. 96, 3.2/.97, 3.5/. 98, 3. 9/.99, 4.6/.995, 5. 3/.998, 6.2/ 
. 999, 7/.9997, 8/ = 
2 FUNCTION RWN2, C29 

0, -3/.012, -2.25/.027, -1. 93/.043,<1. 727. 062, -1.54/. 084, -1.38/ 

Aisa, _1.26/. 131, -1. 12/. 159, -1/. 187, -. 89/. 23, - 94). 267, =. 6e/ 
334, -.43/. 432, -.17/.5,0/.568,. 17/, 666, .43/.732,.62/.77,. 74/ 
.813,.89/. 841, 1/. 988, 2. 2o7i,37 


me NW NH 





. 3 FUNCTION P1,D2 
- 2, 10/3, 18/ 
4 FUNCTION RN3, D2 
. 750, 2/1. 000, 3/ 
- 4 FVARIABLE 16+4*FN2 
- 1 GENERATE 3+?! 
2 ASSIGN 1, FN4 
3 TEST L Q2,2,5 
4 TRANSFER ,7 
5 TABULATE Pl 
6 TERMINATE 
7 QUEUE 2 
8 ENTER 2 
9 DEPART 2 
10 ADVANCE  FN3, FN1 
11 LEAVE 2 
; 12 TRANSFER ,5 
13 GENERATE 180 
14 TERMINATE 1 
: START 1 
END 
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INPUT FILE WITH X-VECTOREQUIVALENT OF SAMPLE PROBLEM 
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X-VECTOR STATUS AT END OF SIMULATION RUN 
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PROGRAM USED TO LOAD X-VECTCR AND CALL SIMULATION ROUTINE 
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ACLOCK 
ALTBLO 
AVAIL 


BASE 


CBLO 
COMVAL 
CQUE 


CREATE 


GSTO 
CTRA 
DELAY 


DTIME 


FTRA 


FREWDT 
GIVEUP 


INST 


APPENDIX B 


DICTIONARY OF PROGRAM VARIABLES 


Absolute Clock Time. 
Alternate Block Number. 
Amount of Storage currently available. 


Value to which spread will be added in GEN- 
ERATE and ADVANCE blocks. 


Pointer to current Block being processed. 
Comparison value used ina SELECT Block. 


Pointer to current Queue being processed. 


Time from current clock time until Transaction 
due to leave FEC. 


Pointer to current Storage being processed. 
Pointer to current Transaction being processed. 
Time delay caused by ADVANCE Block. 


Time differential between clock and last entry 
into a Queue or Storage. 


Stack for evaluating floating point arguments. 


Pointer to first Transaction on list of unused 
Transactions. 


Width of frequency interval ina Table entity. 
Number of Storage units being made available. 


Pointer to. current block element being processed. 
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INTDEC( 


HIGH 


JBL 


KDEX 


KMODE 


LOW 


MAXCTS 
MAXX 


MD FR 


MODE 


MRNG 


NEWBDT 


NEWPRI 


NEXTBDT 


NEXBLO 


NOUPD 


OLDBDT 


OLDPRI 


) 


Indicates whether corresponding OS/FOS value 
is integer or floating point. 


Ending entity number to be used inSELECT 
Block search. 


Number of Block whose routine is being executed. 


Index used to keep count of number of arguments 
specified for a Block. 


Special variable used to determine a unique 
routine to be performed fora given set of 


conditions. 


Beginning entity number to be used in SELECT 
Block search. 


Maximum Contents of Storage or Queue. 
Maximum number of X-Vector elements. 


Modifier used in TEST, SELECT and GATE 
blocks. 


Indicates whether simulation will begin in Clear, 
Reset or Restart mode. 


Number of points ina Function. 


Block departure time of Transaction being 
merged into FEC. 


Priority of Transaction being merged into CEC. 


Block Departure Time of first Transaction on 
FEC. 


Next Block number Transaction is to enter. 


Indicator used to merge a Transaction into CEC 
without updating Clock. 


Block Departure Time of Transactions being . 
checked in FEC. 


Priority of Transactions being checked in CEC. 
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OPCODE 
: os( ) 

PCONTS 

PNTA 


PNTB 
. | PNTF 

PNTS 

PNTT 


PNTX 


SCFLAG 
SE 


SF 





SNE 


SNF 


SRED 


Operation Code indicating type of Block. 
Stack for evaluating integer arguments. 
Present contents of Queue. 

Pointer to last entry on return address stack. 


Pointer to last evaluated argument on OS/ FOS 
stack. 


Pointer to first element of frequency intervals 
for a Table entity. 


Pointer to first element of Slope values of 
current Function entity. 


Pointer to Transaction being compared with 
current Transaction. 


Pointer to first element.of independent Variable 
values of current Function entity. 


Pointer to first element of Function values of 
current Function entity. 


Return Address Stack used to temporarily hold 
location of next arguments to be evaluated. 


Status Change Flag. Set equal to 1 when 
certain entity attributes change and 0 otherwise. 


Pointer to top of push-down chain for empty 
Storage condition. 


Pointer to top of push-down chain for full 
Storage condition. 


Pointer to top of push-down chain for not empty 
Storage condition. 


Pointer to top of Push-down chain for not full 
Storage Condition. 

Pointer to top of Push-down chain for reduced 
contents Storage condition. - 
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SSIND 


“ TASGN 
TBLNO 
TEST 


TPAR 


TRCNT 


TRSIZE 


TYPARG 


TYPX 


~iry 


Scan Status Indicator. True if Transaction active 
in CEC. False if Transaction inactive in CEC. 


The value being assigned in an ASSIGN block. 
Number of Table being tabulated into. 
Value of pointers associated with delay chains. 


Parameter being assigned value in an ASSIGN 
Block. 


Indicates the Transaction number being assigned 
a new Transaction. 


Indicates size of X-Vector printout. 


Indicates whether search argument of a Function 
entity is integer or floating point. 


Indicates whether independent variable of Function 
entity is integer or floating point. 


Indicates whether function values of Function 
entity are integer or floating point. 


Indicates upper limit of lower interval for a 
Table entity. 


Number of units entering or leaving a Queue. 
Current contents of Storage. 
Number of Storage units required. 


Weighting Factor for use in Tabulate block. 
Normally equal to 1. 


Holds statement number for use in returning 
from various routines. 
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Qc 


UE) CALL ERRCR(500+0+6&9999) 


GOTO 190 
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TC 


CALL ERROR(435:0+69999) 


TO 


CALL ERROR (435 10 ,£9999) 
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FN 


FUNCT) CALL ERRGR(507+0s6&9999) 
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as 


) eFOS(PNTB+1)) 


9 
0 
OS(PNTB)*CS(PNTB+1) 
FOS (PNTB)*FOS (PNTB41 ) 
OS (PNTB)/OS(PNTB+1) 
FOS(PNTB)/FOS (PNT B+1) 


MOD(OS(PNTB),OS(PNTB+1)) 
AWOD(FOS(PNTB) »FLOAT(CS(PNTE+1))) 
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